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Abstract: We study a worldline approach to quantum field theories on flat mani- 
folds with boundaries. We consider the concrete case of a scalar field propagating on 
]R + x WL ^ 1 which leads us to study the associated heat kernel through a one dimen- 
sional (worldline) path integral. To calculate the latter we map it onto an auxiliary 
path integral on the full W D using an image charge. The main technical difficulty 
lies in the fact that a smooth potential on R + x M D_1 extends to a potential which 
generically fails to be smooth on M. D . This implies that standard perturbative meth- 
ods fail and must be improved. We propose a method to deal with this situation. 
As a result we recover the known heat kernel coefficients on a flat manifold with 
geodesic boundary, and compute two additional ones, A% and Ar. The calculation 
becomes sensibly harder as the perturbative order increases, and we are able to iden- 
tify the complete Ar_ with the help of a suitable toy model. Our findings show that 
the worldline approach is viable on manifolds with boundaries. Certainly, it would 
be desirable to improve our method of implementing the worldline approach to fur- 
ther simplify the perturbative calculations that arise in the presence of non-smooth 
potentials. 
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1. Introduction 

The worldline approach to quantum field theory (QFT) has been extensively used 
to simplify calculations of many physical quantities like amplitudes, anomalies and 
effective actions, see [0 for a review. Recently this method has been extended to 
include the couplings of particles of spin 0, 1/2 and 1 to external gravity [|, |J 
together with some new applications 0. Extensions to higher spin fields is also under 
study 0. In this paper we investigate the possibility of using a worldline approach to 
quantum field theories defined on manifolds with boundaries. Previous investigations 
of the worldline path integral approach on manifolds with boundaries have already 
been carried out by using numerical Monte Carlo simulations and successfully applied 
to Casimir energy calculations [|7|, ||. 

QFTs on manifolds with boundaries have always attracted a certain interest, with 
applications ranging from Casimir energies to critical phenomena near boundaries, 
including worldsheet approaches to open strings. More recently, they have drawn 
attention in the context of the Horava-Witten theory |J and in the brane world 
scenarios, as for example in [[10] . Thus, it seems useful to investigate and develop new 
calculational tools that may be used to study the properties of QFT with boundaries. 

We investigate here an analytical approach based on first quantization, where a 
path integral over the coordinates of the field quanta is introduced. To reach directly 
the heart of the problem, we consider the case of a scalar field defined on a flat 
space with geodesic boundary, the D-dimensional manifold M. = M + x We 
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use euclidean conventions and consider a generic potential 17(0). The QFT action is 
given by 

S QFT [<f>] = I </"•'■( ^<WrM !'{<>)) (1.1! 



A(^ M 0^0 + C/(0) 



and after specifying suitable boundary conditions (e.g. Dirichlet or Neumann) one 
may compute the corresponding one-loop effective action, which can be represented 
as 



log Det-a [-D + U" 
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The path integral in the last line is evaluated with periodic boundary conditions 
(PBC) and contains the worldline action 



S[x] 



<1t ( + U"(0) 



1.3) 



which produces the differential operator — □ + 17" (0) as quantum hamiltonian. The 
scalar field 0(x) is defined on M. and, as a consequence, the paths x^[t) are restricted 
to lie on J\4, see figure [L| Thus, one has to develop tools to calculate worldline path 
integrals with such a restriction on the path integration variables. 




Figure 1: Quantum paths from the ini- 
tial point A to the final point B, in the 
flat space Ai with boundary dM. 
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Figure 2: The image charge method: 
the bouncing path of figure 1 is mapped 
onto the path A — * 1(B). 



General constructions of path integrals on target spaces with boundaries can 
be found in several textbooks, see for example |TT|, |T2j]. However, the time slicing 



constructions presented there do not seem to assist in practical perturbative calcula- 
tions. One would like to have a reliable scheme to perform perturbative calculations 
in the presence of the boundary. For this purpose we find it convenient to introduce 
an image charge and map the original problem to a new one defined on the full space 
R D , see figure 0. The major difficulty now arises from the fact that the potential 
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extended on the full space may be non-smooth at the fixed points of the Z2 identi- 
fication that reintroduces the boundary. In this paper we develop a scheme to deal 
with this situation, and use it to compute perturbatively the matrix elements of the 
heat kernel (y\e~ T( >~ n+u " ( -^\x) together with its trace, which appears in ( |1.2|) . The 
non-smooth part of the interaction potential on M. D is going to be due to the pres- 
ence of certain step functions. Our scheme implements these step functions, which 
contain the quantum field and sit inside correlation functions, by transferring their 
effect to the limits of integration of suitably inserted completeness relations. 

Several methods for computing the trace of the heat kernel on manifolds with 



boundaries have been studied and used in the literature, as for example in [14 



15 , |i~6j |i~7| , |i~8l , |19 . For a review and additional references one may consult |20fl . The 
method which is closer to the present one is presumably that developed by McAvity 



and Osborn fl7|| , where the DeWitt ansatz |2T| for the heat kernel is generalized 
to the case with boundaries. The novelty of our approach is that we use a path 
integral to calculate directly the heat kernel. An advantage of the path integral 
method is that it is quite intuitive and flexible, allowing to study the complete heat 
kernel and more general correlation functions, though we shall not investigate the 
latter here. Originally we also expected it to be quite an efficient method, but we 
find that this expectation is not fully realized, and presumably it can be improved by 
devising alternative and more efficient ways of calculating the perturbative expansion. 
Nevertheless, in this paper we are able to reach higher perturbative orders for the 
trace of the heat kernel than those we could find in the literature [18, 2(J, namely A3 
and Ar, see eqs. ( |3.27|) and Q3.28| ), but to identify the full A 7 we needed to use the 
result of a suitable toy model. In general, higher orders of the perturbative expansion 
produce multiple integrals that are difficult to compute explicitly, even though they 
could be computed numerically, perhaps using the methods of refs. J?], We use 
the toy model to identify precisely the value of one of these integrals. 

After reviewing in section 2 the path integral approach to the heat kernel on 
flat spaces without boundaries, we discuss in section 3 the case with a geodesic 
boundary, presenting first a brief description of our method and then its application 
to the half line and to R+ x E° _1 . We give the heat kernel coefficients A n with n 
integer and half integer up to n = ~. We use a toy model to identify the value of 
a numerical coefficient entering the complete expression for Ar. Finally we present 
our conclusions and an outlook in section 4. 



2. Heat kernel and path integrals 

Let us start reviewing the standard definition of the heat kernel and the way a path 
integral can be used to compute it. This will also serve to define our notations. 

Given the differential operator corresponding to the hamiltonian of a nonrela- 
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tivistic particle with unit mass 

H x = ~Vl + V(x) (2.1) 

where is the laplacian in the cartesian coordinates x^ of flat space ~R D and V(x) 
an arbitrary potential, the heat kernel K(x,y;f3) is the solution of the Schrodinger 
equation in euclidean time (the heat equation) with a particular boundary condition 
at time (3 = 

-^K(x,y;(3) = H x K(x,y;(3) 

K(x,y;0) = 5 D (x-y) . (2.2) 

The heat kernel is formally given by the matrix element of the evolution operator 
e _/3H , which in the Dirac bracket notation reads 

K(x,y;P) = {y\e- f>H \x) (2.3) 

and corresponds to the transition amplitude for the quantum system to start at 
point x^ and reach point y^ in an euclidean time (3. After setting f3 — it it gives the 
transition amplitude corresponding to the Schrodinger equation in real minkowskian 
time t, but we prefer to use the Wick rotated version with euclidean time 1 . 

Various methods for computing the heat kernel from the Schrodinger equation 
with arbitrary potentials V(x) have been developed in the literature. One of the 
most useful is the path integral approach which produces K(x,y;(3) as a sum over 
all paths x^(t) linking the initial point x^ to the final point y^ in a time (3 



px{P)=y 

K(x,y;f3) = Dx e~ s[x] (2.4) 

Jx(0)=x 



where the euclidean action S[x] in the exponent is the one leading to the hamiltonian 
H x , i.e. 

S[ x ] = J dt fy^xW + V(xj) . (2.5) 

The symbol Dx indicates the path integral measure 

Dx= Yl d D x(t) (2.6) 

0<t</3 

where the range of x*(t) is the full R D for any t. 



1 Substituting — > 2T and V(x) — » \U{<f>{x)) one finds the heat kernel needed for the scalar 
field d> described in the introduction. 
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In general, it is not known how to compute the path integral for arbitrary po- 
tentials V(x). A standard approximation method, useful in many contexts, is to 
compute the path integral in a perturbative expansion for small propagation time 
(3 and small distances £ M = y^ — x M . It is often useful to introduce a rescaled time 
r = t/(3 and write the action as 

S[x] =^jdr fyr&V + P*V(xf) (2.7) 

where the dots now refer to r derivatives. Then one may split the action into a free 
part plus an interacting one 

S[x] = S [x] + S int [x] 

So[x] = ^ [ drU^x" , S int [x) = - [ dr(3 2 V(x) . (2.8) 
P Jo P Jo 

The path integral for the free part Sq[x] is exactly calculable, and the interaction part 
can be treated as a perturbation. The terms generated by the interaction potential 
V(x), assumed to be smooth, can be computed by Wick contracting the quantum 
fields and evaluating the emerging Feynman diagrams. This method of computation 
gives an answer of the form 

K(x, y; (3) = e~ 3 ^ n(rr, y; 0) (2.9) 

(27rp) 2 

where x M (r) = x M + (y M — x fJl )r is the classical trajectory of the free equations of 
motion derived form Sq[x], so that So[x] = j^(x — y) 2 - The term {2irf3)~ 2" is the one 
loop correction which gives the correct normalization of the free path integral, and 
Q(x,y;/3) contains the perturbative corrections due to a nonvanishing potential V. 
The perturbative calculation in terms of Feynman diagrams produces an answer of 
the form 

00 

n{x,y;P)~Y,"n(?,V)P n ( 2 - 10 ) 

n=0 

where ao(x,y) = 1 in flat space. The coefficients a n (x,y) are often called Seeley- 
DeWitt coefficients and correspond to a (n + l)-loop calculation on the worldline 2 . 
Thus we see that the DeWitt ansatz ( |2.9| ) - ( |2.10| ) for solving the heat equation 
emerges naturally from the path integral. 



2 This is generically true in curved space, as propagators go like (3 and vertices like with 
(3 considered as the loop counting parameter. However in flat space the potential V is the only 
source of vertices and contains an extra power of f3 2 which offsets this counting: each vertex with 
V increases the loop counting by 2 because of the factor /3 2 . 
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For example, assuming a smooth potential V, one may compute perturbatively 
to order (3a (counting £ M = y^ — x M ~ yf]3) the transition amplitude 



K(x,y;P) 



{2ti(3) 



x exp 



2(3 2 



, x i + -z-d + — tt-d) 2 + — tt-d) 3 )nv(x 

4! V 2^ 20 v ^ ; 30 vs ; 1 



+ 4[( 1 + ^' 9 )( (Wx))2 



10 



□V(x) + 



(2.11) 



We have given the result in an exponentiated form since it is simpler to compute 
the connected worldline diagrams only Note that in this compact notation one 
considers [£, d] = 0. It is sometimes useful to present the result in a symmetrized 
form as follows 



K(x,y;(3) 



e W 



{271(3) 



exp [/?( - v + l^n^ - ^e^rr^^^^^y 



+ 



5! 



(2.12) 



where we have introduced the notation V = h(V(x) + V(y)), etc., so that the result 
is written in a more compact form and the symmetry x y is manifest. 

Expanding the exponent one reads off the expansion in £ of the Seeley-DeWitt 
coefficients ao, ai, 02, 03 

a (x,y) = 1 

ai(x s y) = -V + ^eCd^V - I^-^^^^F + 0(£ 6 ) 



a 2 (ar,y) 



1 



1 



12 



(2.13) 



a 3 (x, y) = + ^ ((<W 2 + 2V - -n*V) + 0(£ 2 ) . 

Their values at coinciding points y^ = x^ 1 (i.e. ^ = 0) are easily obtained 

ao(x, x) — 1 
ai(x, x) = — V 

a 2 (x,x) = V-ln^ 

As for higher orders, one may find in p2 the heat kernel coefficients a n (x, x) (modulo 
total derivatives) up to n = 8 included, a result which has been obtained using a 
worldline approach. 



(2.14) 
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3. Heat kernel and path integrals on a flat space with bound- 
ary 



3.1 The method 

Now let us consider the case of a space with a boundary. For definiteness we consider 
the flat space Ai = M+ x with coordinates x^ = (y, z l ) where < y < oo and 



z l G M D 1 . This space has a boundary dAi = ~R D 1 located at y = 0, see figure |I[ 
We wish to obtain the heat kernel for the operator H x again as a path integral 

rx(l)=x 2 

K(x 1 ,x 2 ;f3) = Dx e~ s[x] (3.1) 

Jx(0)=x 1 

with the same action S[x) as in (p.7|) , but now summed over all paths x M (r) that lie 
in Ai. In the presence of a boundary one must impose boundary conditions on dAi, 
and for simplicity we assume either Dirichlet boundary conditions, K(x, X21P) = 
for x G dAi, or Neumann boundary conditions, d y K[x, x 2 ; 0) = for x G dAi where 
d y is the derivative normal to the boundary. 

It is evident that on such a manifold the free action Sq has two local minima 
corresponding to the straight classical path X\(t) joining x\ to x 2 and to the path 
x 2 (r) which bounces once at the boundary, see figure ^. The path integral can be 
computed by summing over all fluctuations around these two classical paths. This 
produces an answer of the form 

K(x u x 2 - (3) = (e- So[Sl] Qi(x u x 2 - (3) + 1 e~ So ^ Q 2 (x u x 2 - (3)) (3.2) 

(27T/3) 2 v J 

where the Qi(xi,x 2 ; (3) contain the perturbative loop corrections around the classical 
paths Xj(r), and 7 is a relative phase between these paths which depends on the 
boundary conditions chosen on dM.. This formula reproduces the generalized DeWitt 
ansatz introduced by McAvity and Osborn to study the heat kernel on spaces with 
boundaries [ l7| . In the following we describe a method for computing the path 
integral and obtain perturbatively the functions £li(x\, x 2 ; (3). In particular we will 
focus on the calculations of the trace of the heat kernel 

Tre- fiH =f d D xK(x,x;(3)= [ Dx e~ s[x] (3.3) 

J M JPBC 

where PBC denote periodic boundary conditions for the trajectories on Ai. 

For evaluating the path integral over Ai it is convenient to use the method of 
images and double the space to obtain the full M. D . Trajectories that bounce at the 
boundary are traded for trajectories that continue off the boundary with a reflected 
path, see figure |j. A trajectory extended on the full 1R D (dashed line in figure |4]) 
contributes the same weight in the path integral as the original one (solid line), 
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Figure 3: Classical paths in the flat 
space with boundary A4. On the vertical 
axis we have placed the time r. 
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Figure 4: A bouncing quantum path 
and its image on M> D . 



provided the potential is extended on Mr as an even function under the reflection 

y^-y 

V(x) - V(x) = 6(y)V(y, z) + 6(-y)V(-y, z) (3.4) 

where 9 is the Heaviside (step) function. 

Dirichlet or Neumann boundary conditions are special cases of the generic ansatz 
(|3.2|) : they correspond to setting 7 = — 1 or 7 = +1, respectively. Hence, the 
transition amplitude for such cases may be expressed as 

K M (x 1 ,x 2 ;P) = K r d(x 1 ,x 2 ;P)tK r d(x 1 ,x 2 ;P) . (3.5) 

The heuristic interpretation for 7 = — 1 is quite clear. In fact, notice that the second 
term in ( |3.5| ) involves paths that cross the boundary at least once. Hence, it is easy 
to convince oneself that - due to the symmetry of the potential - the only effect of 
such a term is to cancel all those paths in the first term that hit the boundary. In 
fact, setting y 2 = (all paths now hit the boundary) leads to a vanishing kernel: 
therefore the wave function evolved with the kernel 

ip(x 2 ,P) = / dx 1 K M (xi,x 2 ;P)ip(x 1 ,0) (3.6) 
Jm 

satisfies Dirichlet boundary conditions. For 7 = +1 we know of no simple heuristic 
interpretation; however it is clear that differentiating ( |3.5|) with respect to y 2 and 
setting ?/2 = leads to a vanishing result: the normal derivative of the wave func- 
tion vanishes at the boundary. Hence, 7 = +1 corresponds to Neumann boundary 
conditions. 

The problem has thus been reduced to a path integral computation on a bound- 
aryless manifold. However the path integral on IR D with the potential V is not easily 
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computable by Taylor expanding the potential and using Wick contractions because 
of the step function 6(y) contained in V. Even assuming that V(x) is smooth enough 
to admit a Taylor expansion, V(x) in general is not. 

Nevertheless, as a test, one could try to insert V into ( |2.11| ), which was obtained 
by assuming smoothness of the potential. Evaluating the derivatives contained in 
( [2 on V would then produce delta functions and derivatives thereof, which use 



seems rather unwarranted and difficult to implement consistently. 

Thus one has to study methods to compute perturbatively path integrals that 
contain step functions of the quantum fields. A method to address this problem is the 
following one. Whenever a step function containing the quantum field q(r) appears 
inside correlation function, such as (0(g(ri) — a)), one may split the path integral 
into two parts using a completeness relation at time T\ (namely I = ^ dqi\qi)(qi\) 
to transfer the constraint required by the step function into the limit of integration - 
namely J a °° dq\. Each of the two remaining parts are standard path integrals without 
step functions, but have boundary conditions linked by the J a °° dqi integration. 

Next we will exemplify this method to compute the heat kernel expansion for 
a particle on the half line. This D = 1 case is enough for our purposes, since it 
contains the essential information. At the end of the section we shall reconstruct the 
complete results for generic D. 

3.2 Particle on the half line 

We consider here the heat kernel on the half line R + and denote by x its single coordi- 
nate. For Dirichlet/Neumann boundary conditions the expression for the transition 
amplitude reads (in an obvious notation) 

{x 2 \e~ pH \xi)^ + = (x 2 , P\x u 0) M+ = (x 2 , P\x!, 0) R =F (x 2 , P\x u 0) k (3.7) 

where x 2 = —x 2 is the coordinate of the image charge. We now proceed to compute 
the path integral for the two terms on the right hand side, as explained earlier. 

As usual in path integral computations one splits the quantum paths into a 
classical path x c /(r), satisfying the equations of motion, and quantum fluctuations 
q(r). We associate the boundary conditions with the classical path so that the 
quantum fluctuations g(r) have vanishing boundary conditions 

x(t) = x c i(t) + q(r) = x x + {x 2 - x\)t + q(r) , (3.8) 
q(0) = q(l) = . (3.9) 

We treat the potential V as a perturbation for the free particle. Hence, the path 
integral for the full line, appearing in (|3.7| ), is given by 



(x 2 ,P\x u 0)* = e~ So ^ / Dq e- 5 °Mexp (-(3 dr V(x d (r) + q(r)) ) (3.10) 

Jq(0)=0 V JO 
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pq(l)=0 



+ 



P 



2 /•! 



Dq 1 - / dn V(x d ( Tl ) + q(n)) 

g(0)=0 V JO 



d n / rfr 2 V(x d (ri) + ^(^(ra) + g(r 2 )) + 0(1/ 3 ) 



where So is the action for the free particle, V" = 0. The first term in the expansion 
simply yields the path integral normalization 

1 



?(i)=o 

Dq e~ So ^ = 
9(o)=o (2tt/3)5 

whereas the rest contributes to the perturbative corrections to the heat kernel. 












q(0) q(0) 



—X d (Xi) Q 



Figure 5: Quantum paths constrained Figure 6: Quantum paths constrained 
by 6(xd(Ti) + q(n)). The time r runs by 6(-x d (Ti) - q(ri)). 
along the vertical axis. 

We have managed to compute the short-time perturbative expansion of the tran- 
sition amplitude up to order considering x\ — %2 ~ y/p. In the following we only 
give a detailed description of the contributions coming from a single insertion of the 
perturbation V, since this is enough to exemplify the method discussed previously. 
In principle this method can be applied to all orders in V. 

To compute the aforementioned contributions we first extract the integral over 
the worldline (namely dn) out of the path integral; we will perform it at the end. 
Then, using (|3.4j ), we have 



OMn) + qinWixain) + 5 (n)) 

+6{-x c1 {t 1 ) - g(r 1 ))y(-x d (r 1 ) - g(n))l . (3.11) 



9(1)=0 

Dq e- So[q] 

q(0)=0 



Note that the constraints for the quantum field q{r), associated with the Heaviside 
functions, are localized at time r = T\. Each of the two terms in ( 3.1 1|) can thus be 
manipulated by splitting the path integral into two parts, as follows. In the first term 
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the constraint acts as depicted in figure [5|: we can thus compute the path integral as a 
convolution between a free path integral with boundary conditions q(0) = 0, q(r x ) = 
y and a free path integral with boundary conditions q{r x ) = y, q(l) = 0, where y is 
integrated over and satisfies the constraint y > — £ c «( r i) because of the step function. 
We thus have 

[ qil) Dq e~ s °M 6(x d (n) + q(n)) V{x d {r x ) + q{r x )) 

Jq{0)=0 

/°° ( ri( T i)=y \ ( /•9( 1 )=o 

dy I Dq e"^ 1 V(x cl (r x ) + y) / Dq e~ s ^ 

•*el(n) \Jq(0)=0 I \J q ( T1 )=y 



30 g 2/871(1-7-].)" 

dy Vfofa) + y) . (3.12) 

*d(n) 2n/3y/T 1 {l - T X ) 

In a similar fashion, in the second term of (ft.ll ) the constraint acts as depicted 
in figure ||, and the path integral reads 
?(i)=o 

Dq e- So[q] 6{-x d {n) - q(n)) V{-x cl {n) - q(n)) 

?(0)=0 

~ x cl{Tl) g _ 2/3T 1 (l-T 1 )2 /2 

dy V(-x d ( n ) - y) . (3.13) 

27T/5a/ti(1 - n) 

Hence, the complete contribution of a single V insertion to the whole line path 
integral can be written as 



-3!cl(r) y 2 

+ / dy e" 2 ^ 1 —) 



V(x d (r) + y) - V(-x d {r) - y) , (3.14) 



where the constraint only sits in the integration limit of the second term. Higher 
order contributions in V can analogously be split into one unconstrained part where 
the "quantum field" y = q{r) may run over the whole real axis, and parts that 
depend upon the constraints: the latter are parity-odd. A few observations are in 
order. 

• For an even potential V the parity-odd terms vanish identically. In such a case 
it is possible to write the heat kernel as a power series in integer powers of j3 and 
X2 — X\. In fact, upon Taylor expanding the potential about the initial point x x , the 
unconstrained gaussian integral over y singles out only integer powers of j3, and one 
reproduces the expression Q2.12j ), with x — x x , y = xi and £ = Xi — x x , plus an 
image charge contribution. We can thus write down the final expression for the heat 
kernel on the half line 



e w 



{x 2 ,P\x x ,0) M+ = — — yexp 

(271/5)2 



2 



-w - (1 - e/p) 
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I? 



3 



YT^d'V (1 - 2f//3 + 2£ 4 //? 2 ) + ^dV' + 0(/? 4 ) 



e w 
(2tt/3)2 

2-5!' 



exp 



2 ■ 3! 



0HM l _ 2^7/3 + 



+ ^dv 2 + o(p A ) 



(3.15) 



where £ = — x 2 — X\. The expansion above clearly corresponds to the McAvity- 
Osborn ansatz (|3.2|) , with O^'s both expressed as integer power series in j3. In the 
present case with even potentials, it is clear that the heat kernel is also correctly 
reproduced by the conventional perturbative calculation of the path integral through 
the Wick's theorem, as discussed in section 0. 

• For a generic potential V, the constraint-independent part of the expansion yields 
again the expression (p. 15 ), except that now V = h(V(— x 2 ) + V(xi)) must appear 
in the second exponential in place of V = ~(V(x 2 ) + V(a;i)), and so on. However, 
the constraint-dependent part now is not vanishing and does not seem to be naively 
expressible as a power series in (3 and x 2 — x±, as the latter appears nontrivially in 
the integration limits. This complication corresponds to the difficulty of giving a 
reliable interpretation of the Taylor expansion of the Heaviside function, pointed out 
in the previous section. 

• The method described above allows to compute the perturbative expansion in (3 
of the partition function 



Tr e~ 0A 



dx (x,P\x, 0) R+ 
dx (x, f3\x, 0)k =f / dx {— x, (3\x, 0) 



(3.16) 



for which we need to specialize expression ( |3.14|) , and analogous higher-order contri- 
butions, to 



W x i 
and 

(ii) xt 



%2 



x 



Xcl(T) 



X 



-x 2 



X 



x(1-2t) 



Thus we can proceed to compute the small time expansion of the partition function 
up to order (3^\ 

(i) In this case, the constraint-independent parts of the potential insertions, such as 



the first term of ( 3.14 ), yield bulk contributions, whereas the constraint-dependent 
parts, after Taylor expanding the potential about the boundary, single out boundary 
terms. Hence, after some tedious algebra, we have 



dx (x, P\x, 0)j 



(27r/3)i 



dx exp( -PV(x) - f3 2 2 ^ J + (3 3 K A >1 
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2-5! 



P 2-3! 



V 2-3! 2-5! J 



+ 



(3.17) 



(ii) For such a case, the presence of the non trivial classical action So[a; c j] = 2x 2 / (3 
allows for a Taylor expansion of the potential about the boundary location: the 
integral over the manifolds is convergent and gives 



OO 

dx (-x,/3\x,0) R = - — 

(271-/3): 




1 - 07(0) 



,3 / ^ 3 (0) , 3(<97(0)) 2 _ <9 4 V(0) \ 
P V 3! 2 4 2 7 y 



^7(0) 



+ /? 3 



/ 797(0) _ 9¥(0)\ 
V 4 2-4! J 



+ 0((3 4 



(3.18) 



Combining ( |3.17| ) and ( |3 . 1 8| ) we get the final result 



Tr e 



-/3H 



dx (x, P\x, 0)i 



(27T/3) 



J dx exp^—/3V(x) — (3 



3 d 4 V(x)\ f^ Q 7(dV(0)f 



2 d 2 V(x) , a3 (dV(x)) 2 



2 ■ 3! 



+ /3 J 



2-5! 



+ 



'tt/3 



4! 



, & (_YM , ^ 2 7(Q) , (^(Q)) 2 _ 5Mo)^ 

+P I 3! 2 3 4! 2 7 J 



-11 797(0) f 2 | <9 3 7(0) \ 
2 J 3! + \-3j 5! J 



(3.19) 



0(/3 4 ) 



This is the expansion of the partition function with Dirichlet /Neumann boundary 
conditions, valid up to order /?2, where the upper coefficient in the symbol {:} refers to 
Dirichlet boundary conditions, and the lower one to Neumann boundary conditions. 
At this point the integrated heat kernel coefficients for the half line can be read off 
immediately. However, it is not difficult to reconstruct such coefficients for arbitrary 
D. Defining 



Tr e 



(3.20) 



13 



we obtain 



A, 



Ai 



= / 

J A 



A 1 



2 



A, 



As 

2 



M 



^ § Jqm 



Jm 



V 



M 



JdM 



V 2 uv 



12 



+ 



1 

-2 



dM 



dyV 

3! 



'7T 



92 a? 



2 1 8 12 1 



+ 



dM 



V 3 



4! 



+ 



vnv u 2 v 



-1 1 VdyV 
3! 



+ 



2 • 3! 
2 



2-5! 



At 

2 



7T 

8 ./aM 



y3 
3! 



-3 
4! 



5! 



3 

-7 



2-5! 



<9 2 d 2 ' 
1 2 3 12 1 



9, 4 d 2 d 2 
+ 



+ 



a? 



2 7 4-4! 10 • 4! 



+ 



-51 (W 
7 | 2 6 



(3.21) 
(3.22) 
(3.23) 
(3.24) 
(3.25) 
(3.26) 



(3.27) 



(3.28) 



Before concluding this subsection let us mention some caveats of our method. We 
have found it increasingly harder to calculate explicitly the numerical value of some 
terms. On the other hand our method allows to write all such numerical values 
in terms of multiple integrals of finite functions over compact domains which could 
eventually be calculated numerically. Since in the present manuscript we only en- 
countered difficulties for one numerical coefficient in Az, and precisely the one related 
to (d y V) 2 , we decided to fix it by using a toy model, which is discussed next. 



3.3 Toy model: V(x) =ax 

The asymptotic expansion of e~^ H for small f3 can also be obtained from the large- 
A asymptotic expansion of the resolvent (H + A)^ 1 . Indeed, since the trace of the 
resolvent is the Laplace transform of the heat kernel trace, it admits the following 
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representation for large A [20 



Tr (H + A) -1 ~ ^r^ + ^r-^ 



(3.29) 



neN/2 



We use such an expansion to compute the (d y V) 2 term in the coefficient Ar, con- 
sidering a toy model that can be explicitly solved in the half line, namely the one- 
dimensional Schrodinger operator H = — ^d 2 + ax with x G 1R + and a a real constant. 
We can compute the contribution to proportional to (d y V) 2 from the order A~ 4 
term of the asymptotic expansion, since the complete A r for the linear potential 
reduces to 



Ar 

2 



8 2 6 



(3.30) 



with (<9 y l/) 2 = a 2 for a constant k which has to be fixed. 

From the eigenf unctions of H one obtains a simple expression for the kernel 
G(x, x', A) of the resolvent (H + A) -1 



G(x, x', A) = -— {9(x' - x)L(x)R{x') + 9(x - x')L(x')R{x)} 
where R(x) can be written in terms of the modified Bessel function K v (z) 

13/2- 



R{x) = Vax + X Ki/3 ( — — 3q J 



and L{x) is given by 



- ^1/3 I — 3^ — J ^1/3 I 3^ 



for Dirichlet boundary conditions at x = 0, and 



L(z) = v / axTA{v / 2A 



1/3 I 3a 



(2A) 3 / 2 \ T ( [2(arE+A)l 3 / 2 
J V3 I 3~ a 



(3.31) 



(3.32) 



(3.33) 



j 2_ 



A/3 ) A V3 ^ 

f (2A) 3 / 2 \ r 



[2(a:r+A)] 3 / 2 
3a 



+ 



[2{ax+X)f/ 2 
3a 



(3.34) 



-hi, 



(2\f/ 2 \ v _ f[2(ax+A)] 3 / 2 



3a 



#1/3 



3a 



for Neumann boundary conditions at x — 0. is the Wronskian between the 
corresponding L(x) and R(x). 
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From these expressions one can obtain the bulk and boundary contributions to 
the asymptotic expansion of the resolvent trace 

POO 

Tr (H + \y 1 = dxG(x,x,X). (3.35) 
Jo 

From the asymptotic behavior of the modified Bessel functions for large values of 
their arguments, it is simple to see that the boundary contributions of order A -4 and 
proportional to a 2 are given by — A -4 for Dirichlet boundary conditions, and 
j^o? A -4 for Neumann boundary conditions. Comparing these values with ( |3.29| ) 
and ( |3.30| ), we obtain k = —5 for Dirichlet boundary conditions, and k = 7 for 
Neumann boundary conditions. These are the explicit numerical values of some 
multiple integrals, produced by our general method, that we were not able to evaluate 
directly. 

4. Conclusions 

We have developed a path integral method to generalize the worldline formalism to 
manifolds with boundaries. We considered the flat manifold R + x R D_1 and made use 
of an image charge to extend the problem to ~R D . Smooth potentials on M + x M D_1 
extend to potentials on M D which are generically not smooth at y = 0. Our proposal 
to deal with this situation is in principle applicable to any perturbative order, but the 
explicit calculations become increasingly harder as the perturbative order increases. 
It could be useful to develop simpler strategies to evaluate the path integral in the 
presence of non-smooth potentials. This would make it easier to extend the method 
to curved manifolds with boundaries. Nevertheless we have been able to show the 
feasibility of the worldline formalism in the presence of boundaries, and used it to 

calculate two new coefficients, A3 and At, for a flat space hamiltonian with a generic 

2 

scalar potential. On top of using curved manifolds, one could try to extend the 
worldline formalism on spaces with boundaries to include fields with spin and use 
more generic boundary conditions. 
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